Electron and ion thermal forces in complex (dusty) plasmas 
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Expressions for the ion and electron thermal forces acting on a charged grain, suspended in a 
weakly ionized plasma subject to temperature gradients, are derived. These forces can significantly 
affect grain transport in many realistic situations, ranging from complex plasmas in laboratory and 
microgravity experiments, to charged dust in planetary atmospheres and fusion devices. 
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The momentum exchange between different species 
plays an exceptionally important role in complex (dusty) 
plasmas - multi-component systems consisting of charged 
micron-sized ("dust") grains embedded in a plasma [lj- 
yj]. In a weakly ionized plasma, the momentum transfer 
in collisions with the neutral gas "cool down" the sys- 
tem, in particular grains and ions, and introduce some 
damping. In addition, ion and electron flows are usu- 
ally present in plasmas due to global electrical fields (e.g. 
ambipolar or sheath fields in gas discharges). The mo- 
mentum transfer from flowing ions and electrons can be- 
come a dominant factor affecting and regulating static 
and dynamical properties of the grain component. A re- 
markable example wherein the momentum exchange in 
ion-grain collisions plays a crucial role is the formation 
of a void in the central part of an rf discharge in ex- 
periments under microgravity conditions [H-Q. ^ n ^ ms 
case, the momentum transfer rate from the ions stream- 
ing outwards exceeds the electrical force (pointing to the 
center), leaving a void. 

Normally, the momentum transfer is associated with 
the relative motion between the grain and other plasma 
components. The related forces are then called "drag 
forces" (terms friction and resistance force are also rele- 
vant, especially when the grain moves with respect to the 
stationary plasma background). Usually the mean free 
paths of all plasma species (neutrals, ions, and electrons) 
are very long compared to the grain size [~ O(fiia)], 
which corresponds to the limit of large Knudsen number. 
The neutral drag force in this regime was calculated al- 
most a century ago by Epstein Q and is often referred to 
as Epstein drag. The ion drag force is considerably more 
complicated issue. First, in addition to direct ion-grain 
collisions, elastic ion scattering in the electrical potential 
around the grain has to be taken into account. Second, 
the ion-grain interaction under typical experimental con- 
ditions is not weak, which implies that the conventional 
theory of Coulomb scatterin g is inapplicable to describe 
elastic ion-grain collisions |9l llCf. This explains why the 
problem of the accurate determination of the ion drag 
force in complex plasmas has received considerable at- 
tention in the last decade, including various theoretical 
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approaches [91414J , numerical simulations |15Wl7| and ex- 
periments [18l - l23| . Moreover, theoretical results on the 
ion-grain scattering in the regime of strong (Yukawa) in- 
teraction found their application not only in the field of 
complex plasmas, but also in the context of quark-gluon 
plasma [24| , scattering of dark matter particles [25| , and 
conventional strongly coupled plasmas [2f|. Regarding 
the momentum transfer in electron-grain collisions, it of- 
ten plays only a minor role due to the small electron 
mass, although this is not always the case [27j . 

The relative motion between the grain and other 
plasma species is not the only mechanism which can be 
responsible for the momentum transfer. The latter can 
occur even in the absence of net flows of plasma species 
relative to the grain. A relevant example is the ther- 
mophoretic force (terms "thermal force" and "radiome- 
ter force" are also employed in the literature (28|j). It 
describes the phenomenon wherein small particles, sus- 
pended in a gas where a temperature gradient exists 
(but macroscopic flows are absent), experience a force in 
the direction opposite to that of the gradient. Elemen- 
tary consideration of this phenomenon has been already 
given by Einstein (29j and then many others investigated 
this topic in detail (see e.g. Refs. |30l - [33 ] and references 
therein). In the context of complex plasmas, Jellum et 
al. (33l | were apparently the first who recognized the pos- 
sibility to manipulate the particles in gas discharges using 
the thcrmophorctic force. Since then, applying the verti- 
cal temperature gradients to compensate for the particle 
gravity has become a standard technique for controlled 
particle manipulation in laboratory experiments (341437! . 

In complex plasmas, the thcrmophorctic force has its 
counterparts associated with the charged electron and 
ion components, provided the corresponding temperature 
gradients are present. We will refer to these forces as 
the electron and ion thermal forces to distinguish them 
from the conventional drag forces. To the best of our 
knowledge, these forces have not yet been considered in 
the literature in any reasonable detail. The main purpose 
of this Letter is to fill this gap. In the following we derive 
the expressions for the ion and electron thermal forces 
acting on a highly charged grain immersed in a weakly 
ionized plasma subject to temperature gradients. Wc 
analyze the directions and the magnitudes of these forces 
and show that they exhibit interesting and non-trivial 
behavior. Comparison with other forces indicates that 
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the thermal forces can provide important contributions 
to the net force balance in complex plasmas. 

The general expression for the force associated with the 
momentum transfer from a light species a to the massive 
grain at rest is 



m a / vva a (v)f a (y)d 3 v, 



(1) 



where m a , / Q (v) and cr a (v) are the corresponding mass, 
velocity distribution function, and (velocity dependent) 
momentum transfer cross section (a — n,i,e for neu- 
trals, ions, and electrons, respectively). The net momen- 
tum transfer occurs when the velocity distribution has 
some asymmetry (e.g. relative motion). Assuming weak 
asymmetry, we write / a (v) ~ f a o{v) + / a i(v), where the 
symmetric component f a Q is taken to be Maxwcllian 

faO = n Q (m Q /27rT Q ) 3/2 exp(-m Q w 2 /2T Q ). 

Here n a and T a are the density and temperature (in en- 
ergy units) of the species a. The asymmetric component 
f a i, which gives contribution to the integral in ([TJ, de- 
pends on the nature of the anisotropy. In the case of 
subthcrmal drifts with relative velocity u Q {u a < vr a ) it 
reduces to f a i ~ f a0 (vu a /v^ a ), where v Ta = ^/T a /m a 
is the thermal velocity. This ansatz corresponds to the 
conventional calculation of the neutral, ion and electron 
drag forces and has been thoroughly investigated ear- 
lier [2, US, II El H3- If the drift is not subthermal, non- 
linearized shifted Maxwellian distribution function [l4j] 
or another distribution function, appropriate for the sit- 
uation considered [HI, HU, Hsl has to be used in Eq. (QJ . 

The focus of this Letter is on a complementary situa- 
tion when relative drifts are absent, but the momentum 
transfer do occur due to the net momentum flux caused 
by the temperature gradients. In this case, the asym- 
metric part of the velocity distribution function of the 
component a can be approximated as 



fai 



m aK a faO 

n a T 2 



1 



5T a 



vVT Q , 



(2) 



where n a is the thermal conductivity of the species a. 
This form ensures that the self-consistent density gradi- 
ent and/or electric field (for charged species), build up 
in response to the temperature gradient, result in no net 
flux: ja = J vfaid 3 v = (i.e., u Q = 0). In addition, Eq. 
([2]) obviously satisfies the Fourier's law for heat transfer: 
q a = J (m a v 2 /2)vf a id 3 v = — k q VT q . Equation @ is an 
approximation, which is exact only for the special case of 
cx r -4 interactions (for the rigorous mathematical treat- 
ment of non- uniform gases see Ref. 03)- Its accuracy is, 
however, more than acceptable for our present purposes, 
especially in view of further simplifications involved in 
treating electron- and ion-grain collisions. 

Substituting Eq. @ into JT]) we get after integration 
over the angles in spherical coordinates 



Ft a 



16k q VT q 



15v27twt 



x 2 (§ — x) exp(— x)a Q (x)dx, (3) 



where the remaining integration is over the reduced ve- 
locity x = v 2 /2v^ a . Let us first apply this ansatz to cal- 
culate the thermal force acting on a non-charged particle 
in a non-ionized atomic gas - thermophoretic force. In 
this case the momentum transfer cross section is velocity- 
independent, <r n (x) = ira 2 , and the integration yields 
Ft,i = — (8v / 27r/15)(K n a 2 /t'T„)Vr„, where a is the parti- 
cle radius. This coincides with the celebrated expression 
by Waldmann [3(j • The force pushes the particles into 
the region with lower gas temperature. Physically, this 
is because hotter atoms transfer more momentum to the 
grain than the colder ones. 

It is convenient to write the generic expression for the 
thermal forces in complex plasmas as 



(4) 



8V2^ n a a 2 VT a ^ 
F Ta = — 
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For the neutral component (thermophoresis) = 1. 
The factors ^i( e ) account for the electrical interactions 
between the ions (electrons) and the charged grain. They 
depend on the shape of the plasma electrical potential 
around the grain. We use the conventional Debye-Hiickcl 
(Yukawa) form, <j)(r) ~ {Q/r) exp(— r/A), where Q is the 
grasin charge and A is the plasma screening length. In 
low-temperature gas discharges the dominant charging 
process is the continuous absorption of electrons and ions 
on the grain surface. In this case the charge is negative 
and proportional to the product of the grain radius and 
the electron temperature, viz. Q = —z(aT e /e), where z 
is the coefficient of order unity, which depends on plasma 
parameter regime |4ll |. The screening comes from redis- 
tribution of plasma electrons and ions in the vicinity of 
the charged grain. In the absence of substantial ion flows 
and strong nonlinearities in ion-grain interaction (4lT[43l | , 
the screening length is approximately given by the ion 
Debye radius A ~ Ad, = ^/Ti/4ire 2 rii, provided T e ^> T,;. 
Although it is well recognized that the long range asymp- 
tote of the electrical potential can be modified (by e.g . 
continuous plasma absorption on the grain surface [44j — 
or plasma ionization/recombination effects (ill), this 
is expected to affect merely grain-grain interactions, but 
not the momentum transfer from the ions and electrons. 

Electrons. Normally, the grain radius in complex plas- 
mas is much smaller than the plasma screening length. 
In this case, the electron-grain interaction can be called 
"weak" , in the sense that its range - the Coulomb radius 
Rr\ ~ za - is much smaller than the screening length 
A [3| . This implies that the Coulomb scattering theory is 
appropriate to describe electron-grain elastic collisions. 
The general momentum transfer cross section for scat- 
tering in the Coulomb potential Q/r is 



o* a (v) = 4nR 2 a (v)ln 



1/2 



(5) 



where R a {v) = (\Q\e/m a v 2 ) and a = e,i. The maxi- 
mum impact parameter p max is necessary to avoid the 
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FIG. 1: The transition curve separating regions of positive 
and negative electron thermal force in the plane of reduced 
grain size a/X and charge z. 



logarithmic divergence of the cross section. In the stan- 
dard Coulomb scattering theory p max = A, and this 
choice is appropriate for the weak electron-grain inter- 
action 0, ul\ . The minimum impact parameter p m i n is 
zero in the standard Coulomb scattering theory. In the 
considered case, however, some electrons are falling onto 
the grain instead of being elastically scattered in its elec- 
trical potential. Consequently, p mm should be set equal 
to the maximum impact parameter corresponding to elec- 
tron collection by the grasin, p c . The orbital motion lim- 
ited (OML) theory yields (4j| p c (v) = o.y/1 - 2R e (v)/a 
for 2R e {v) < a and p c = otherwise (sufficiently slow 
electrons cannot be collected even in head-on collisions 
with the grain due to electrical repulsion). The momen- 
tum transfer cross section for electron collection is sim- 
ply cf c e (v) = np^(v). Expressing now velocity in terms 
of x = v 2 /2v\^, substituting a e (x) = (J s e (x) + a%(x) into 
Eq. ([3]) and performing the integration with the appro- 
priate limits yields: 



(l + |z + z 2 )cxp(-z)- 
where A e is the electron-grain Coulomb logarithm 



(6) 



A, 



h(x) In (1 



4\ A x< 



dx-2 / h(x)\n(^-l)dx 



(7) 

The function h(x) is defined as h(x) = (| — x) exp(— x). 
Equations (UJ), © and (|7|) constitute tfie expression for 
the electron thermal force. 

For an uncharged particle (z = 0) we have "I> e = 1, as 
expected. In this case the electron thermal force pushes 
the grains in the direction of lower electron tempera- 
ture, similarly to the thermophoretic force. According to 
Eq. ([5]) the contributions from collection and scattering 
are directed oppositely to each other. For a sufficiently 
high charge, the scattering part becomes dominant and 
the force reverses direction. The physical reason is that 
the (Coulomb) scattering momentum transfer cross sec- 
tion quickly decreases with velocity, so that the cold elec- 
trons are more effective in transferring their momentum 
upon scattering. In this regime ($ e < 0), the thermal 
force acts in the direction of higher electron temperature. 



In Figurc[T]thc curve separating the positive and negative 
values of $ e is plotted in the plane {a/X, z). For most 
experimental conditions we expect $ e < 0, although the 
transition line does not seem unreachable [50(. 

Ions. The Coulomb radius of ion-grain interaction 
R l c ~ zar is not necessary small compared to the plasma 
screening length due to the presence of a (normally) 
large factor r = T e /Ti - electron-to-ion temperature ra- 
tio Q . The interaction range can exceed A and consider- 
able amount of momentum transfer can occur for impact 
parameters beyond A. This implies that the standard 
Coulomb scattering approach is inappropriate 0, [TJJ] . It 
makes sense to consider two regimes of ion scattering sep- 
arately. In the regime of moderate ion-grain interaction, 
an extension of the standard Coulomb scattering theory 
is possible by taking into account the momentum transfer 
from the ions that approach the grain closer than A [TJ]. 
This results in p max = A-^l + 2i?,;(u)/A. This approxi- 
mation demonstrates good accuracy for Ri(v) < 5A and 
reduces to p max = A in the limit of weak ion-grain interac- 
tion [Ri(v) <C A]. The impact parameter corresponding 
to the ion collection is p c (v) = ayl + 2Ri(v)/a and the 
corresponding collection cross section is a c (v) = irp 2 (v). 
Combining the contributions from collection and scatter- 
ing yields: 



*i = 1 - jzr-zVAi, (8) 
where Aj is the (modified) ion-grain Coulomb logarithm 

~2x(\/a) + zt 



A, 



h(x) In 



2x + ZT 



dx. 



(9) 



In typical complex plasmas with A 3> a, r 3> 1, and 
z ~ 1, the Coulomb logarithm can be roughly estimated 
as Aj ~ ln(l + P^ 1 ), where /3t = PivrJ = {a/X)zT and 
j3{v) = Ri(v)/X is the ion scattering parameter Q. The 
approach is reliable up to /3t ^5 5. Since the product zt 
is normally quite large, zt ~ O(10 2 ), it follows that (i) 
the scattering provides dominant contribution to the mo- 
mentum transfer; (ii) < 0, i.e. the grains are pushed 
into the region with higher ion temperature. The physi- 
cal reason is again fast decrease of the scattering momen- 
tum transfer cross section with the ion velocity, so that 
cold ions transfer more momentum to the grain. This is 
not a unique example when the force acting on a charged 
grain is directed oppositely to the net ion momentum 
flux. Another example is related to the sign reversal of 
the ion drag force acting on an absorbi ng p article in the 
highly collisional (continuum) limit [4l|, l5ll - [54j , although 
the detailed physics is different. 

In the regime of very strong ion-grain interaction, the 
scattering is characterized by the formation of a poten- 
tial barrier for ions with impact parameters above the 
critical one p„, which considerably exceeds A. For p < p* 
scattering with large angles occurs, which gives the ma- 
jor contribution to the momentum transfer. Relative im- 
portance of momentum transfer from distant collisions 
with p > p„ decreases rapidly with the increase in ion- 
grain interaction strength. The detailed consideration of 



4 



the momentum transfer in this regime can be found in 

Refs. P3, EH. In the limit M v ) > A ( Le - P( v ) > 1 H) 
the total momentum transfer cross section (collection and 
scattering) can be approximated as <7s(u) ~ wp%(v), with 
pi ~ A 2 {In 2 \p(v)\ + 2\n[/3(v)}}. The integration yields 



$ 4 



x 2 h(x) [in 2 (^) 



2 In 



dx. 



(10) 

The upper limit of integration can be chosen as = 
/3t/2 to avoid unphysical regime of negative cross section 
in this approximation. However, due to the presence of 
exponentially small term in h(x), taking x* = oo will not 
produce big errors for /3t 3> 1. Note that the sign of $j 
is changed from negative to positive upon increasing /3t 
(in the considered approximation this happens at /3t — 
120; the exact value is rather sensitive to the functional 
dependence of the cross section on the ion velocity and, 
thus, is subject to significant uncertainty). Physically, 
the sign reversal occurs because scattering in the Yukawa 
potential in the limit of strong interaction tends to that 
on a hard sphere with the radius p*, which only weakly 
(logarithmically) depends on the ion velocity. 

Having derived the expressions for the ion and elec- 
tron thermal forces, let us discuss their role in practi- 
cal situations. As a first example, consider creating a 
gradient of the neutral gas temperature to compensate 
for grain gravity in low ionized laboratory gas discharges 
(thermophoretic levitation) [H, H3|- In this case, the ion 
temperature is likely coupled to the neutral gas temper- 
ature, VXi ~ VT„. The neutral and ion thermal force 
are directed in the opposite directions (in the regime of 
weak and moderate ion-grain interaction). Their ratio 
is \F Tl /F Tn \ — (n l /n n )(a nn /a ln )z 2 T 2 A i , where we take 
into account that n a ~ (n a /n n )(vT a /c an ) in weakly ion- 
ized plasmas (a nn and cr^ are the transport cross sec- 
tions for neutral-neutral and ion-neutral collisions, re- 
spectively) . Assuming Uj n ~ a nn [56[ and A{ ~lwe find 
that F^i dominates over F^ n for n,/n„ > {zt)~ 2 ~ 10~ 4 . 
For a more typical (in laboratory gas discharges) ioniza- 
tion fraction rii/n n ~ 10~ 6 , the ion thermal force is only 
slightly reducing the effect of thermophoretic force. 

The electron thermal force is expected to play more im- 
portant role. Spatially resolved probe measurements of 
the electron energy distribution function in various gas 
discharges have previously shown that significant elec- 
tron temperature gradients can exist [13, He|. For devices 
used in complex plasma research, the results from probe 
measurements [5!| , optical emission spectroscopy [60| , as 



well as from numerical modeling [61], |62[ all revealed gra- 
dients in T e of the order of C(eV/cm). The spatial dis- 
tribution of T e generally depends on discharge geome- 
try, plasma parameters, and can be affected by the pres- 
ence of grains [6l[. Let us therefore assume VT e = 1 
eV/cm and restrict ourselves to the comparison of the 
absolute magnitudes of the electron thermal force and 
other forces acting on a small grain in the bulk of a 
gas discharge. We take the plasma parameter set from 
Ref. [3] used to estimate relative importance of the 
electrical and ion drag forces: argon gas at a pressure 
p = 10 Pa (n„ ~ 2 x 10 15 cm~ 3 ), n e = m = 3 X 10 9 
cm" 3 , T e = 1 eV, = T n = 0.03 eV, a = 1 fim 
(a/A ~ 0.04), and z ~ 3 (estimated from the collisionlcss 
OML theory). From Eqs. © and (J7J we get $ e ~ -22. 
The electron thermal conductivity in a weakly ionized 
plasma is n e ~ (bn e T e /2m e u en ) with v en ~ n n a en VT e 
{&en ~ 10~ 16 cm 2 for T e ~ 1 eV in argon). This results 
in Fre — 2 x 10~ 8 dyne. This force is more than three 
times larger than the force of gravity, experienced by the 
grain of this size (and material density 1.5 g/cm 3 ) in 
ground-based experiments. The equivalent electric field 
for which F Te ~ \Q\E*, is ~ 5 V/cm. Such a field 
would produce a significant plasma anisotropy, character- 



ized by superthermal ion flows (Fig. 4b from Ref. [141]). 
Finally, this magnitude is comparable to the maximum 
value of the ion drag force the grain can experience in 
subsonic ion flows for this set of parameters (~ 5 x 10~ 8 
dyne according to Fig. 4a from Ref. PH)- 

Observation of large grains, trapped in standing stria- 
tions of a stratified dc glow discharge [HI], provides an- 
other example where the electron thermal force can play 
a significant role. The electron temperature is known to 
increase considerably in the head of the striation, which 
can explain that even very massive grains can be con- 
fined there [63[ . Other situations wherein plasma thermal 
forces can play significant role include charged aerosols, 
dust in planetary (e.g., Earth) atmospheres and in fusion 
devices. In the latter case, however, the regime of fully 
ionized magnetized plasma, considered recently [64|, is 
apparently more relevant. 

Altogether, there is a certain confidence that the ther- 
mal forces (especially from electrons) can have significant 
impact on the particle transport in many realistic situ- 
ations. This finding should be properly addressed when 
developing new (and updating existing) numerical codes 
to model particle transport in various plasma environ- 
ments. Expressions derived in this Letter can serve as a 
theoretical basis for further detailed studies. 
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